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ABSTRACT 

Context. The prompt emission of gamma-ray bursts is mostly thought to be produced in internal shocks of relativistic shells emitted 
by the progenitor at different times, whereas the late multi-band afterglow is interpreted as the synchrotron emission of electrons 
swept up by the fireball expanding through the surrounding interstellar medium. The short timescale variability observed in flares 
superimposed on the X-ray/optical afterglow of several bursts recently made possible by Swift, has been interpreted as evidence for 
prolonged activity of the inner engine through internal shocks. Yet, it is not clear whether this applies to all the observed bursts and, 
in particular, whether the bursts exhibiting single y-ray pulses with no short timescale variability at late times could also be entirely 
interpreted as external shocks. 

Aims. We present prompt y-ray, early NIR/optical, late optical and X-ray observations of the peculiar GRB 070311 discovered by 
INTEGRAL, in order to gain clues on the mechanisms responsible for the prompt y-ray pulse as well as for the early and late multi- 
band afterglow of GRB 0703 1 1 . 

Methods. We fitted with empirical functions the gamma-ray and optical light curves and scaled the result to the late time X-rays. 
Results. The //-band light curve taken by REM shows two pulses peaking 80 and 140 s after the peak of the y-ray burst and possibly 
accompanied by a faint y-ray tail. Remarkably, the late optical and X-ray afterglow underwent a major rebrightening between 3 x 10 4 
and 2x 10 5 s after the burst with an X-ray fluence comparable with that of the prompt emission extrapolated in the same band. Notably, 
the time profile of the late rebrightening can be described as the combination of a time-rescaled version of the prompt y-ray pulse and 
an underlying power law. 

Conclusions. This result supports a common origin for both prompt and late X-ray/optical afterglow rebrightening of GRB 0703 1 1 
within the external shock scenario. The main fireball would be responsible for the prompt emission, while a second shell would 
produce the rebrightening when impacting the leading blastwave in a refreshed shock. 

Key words, gamma rays: bursts; X-rays: individual (GRB 07031 1) 



1. Introduction 

The bewildering variety of the long y-ray prompt emission 
profiles of gamma-ray bursts (GRBs) has been puzzling as- 
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tronomers since their discovery (e.g. see Fishman & Meegan 
1995 for a review). Among the most important open issues, two 
are still nowadays debated: first, how long does the central en- 
gine remain active? Does the emission consist of a single episode 
or a temporal sequence of events with interspersed quiescent pe- 
riods? Secondly, after the discovery of the delayed multi-band 
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long-lasting emission called "afterglo w", made possible for the 
first time by BeppoSAX ten years ago dCosta et al.Lll997h . a gen- 
eral consensus on the emission mechanisms of both the prompt 
and the afterglow emission is still missing, even though for the 
latter synchrotron radiation by a population of shocked electrons 
proved to be successful in accounting for a number of spec- 
tral and temporal evolution properties (e.g., Meszaros 2006). 
Gaining clues on these issues may help to shed light on the na- 
ture of the progenitors as well as on the circumburst environ- 
ment. 

The advent of Swift dGehrels et al.L l2004t) has been allow- 
ing the study of multi-band afterglows as early as ~ 10 2 s af- 
ter the prompt event, exploring a previously unknown time do- 
main. Both optical and X-ray light curves have shown unex- 
pected behaviours that still lack a solid agreed theoretical ex- 
planation. Among the impressive discoveries by Swift concern- 
ing the long GRBs, here we mention the canonical behaviour of 
steep-shallow-normal decay characterising most of the early X- 
ray afterglows as well as the occurrence of X-ray flares in ~ 50% 
of them (see Zhang 2007, for a review up to date). 

A number of different interpretations of the shallow decay 
phase experienced by the X-ray afterglows of a sizable frac- 
tion of Swift bursts dTadiaferri et aUl2005l: IO'Brien et al.Ll2006t 
iNousek et all 120061 : IZhang et all 120061) have been proposed. 
They can be broadly classified according to which mechanism 
is invoked: internal (IS; e.g., Rees & Meszaros 1994; Kobayashi 
et al. 1997) or external shocks (ES; e.g., Rees & Meszaros 1992; 
Shaviv & Dar 1995; Fenimore et al. 1996; Dermer & Mitman 
1999). For instance, within the IS model the shallow decay phase 
could be the result of low velocity contrast di stribution wind of 
slow shells emitted soon after the fast ones dGranot & K umar. 

2006) . According to alternative ES interpretations, it could be 
the result of delayed energy injection into the fireball, either 
in the form of freshly ejected material from late engine activ- 
ity ("refreshed shock scenario"; Rees & Meszaros 1998), or by 
the collision of a w ind of low velocity contr ast shells catching 
up with the fireball (Zhang & Meszaros, 2001). More generally, 
the shallow decay could be the result of a "late prompt" emis- 
sion, i.e. the same mechanism responsible for the y-ray pulses 
of the prompt would be at w ork through later shells with de- 
creasing bulk Lorentz factors (Ghisell ini et al.L 120071) . Overall, 
a consensus is still missing, also because of the large variety of 
behaviours observed: in some cases, there is no spectral evolu- 
tion across the break in the light curve marking the end of the 
shallow phase, compatible with the expectations from an exter- 
nal o rigin of the shock s, while for other GRBs the opposite is 
true dLiang et aUl2007l) . 

Likewise, X-ray flares have been interpreted as the result 
of late internal dissipation rather than due to external shocks, 
one of the main arg uments being their short timescales also at 
late times, At/ t «i 1 dBurrows et all 120051; iFalcone et all [2 006; 
Romano etail 120061; IChincarini et all 120071: lLazzati & Pernal 

2007) , whereas an external shock ori gin would require increas- 



ing tim escales (Fenimor e et all[T996l) . although see also lDermeri 
d2007bl) . 

Another debated topic concerns the presence of optical 
brightening contemporaneous with the prompt emission or oc- 
curring in the first 1-2 hours. A variety of mechanisms have been 
proposed to explain it, such as reverberation of the prompt emis- 
sion radiation (GRB 050820A: Vestrand et al. 2006), reverse 
shock (GRB 990123: Akerlof et al. 1999), refreshed shocks 
and/or energy injection (e.g. GRB 021004: Fox et al. 2003; 
GRB 050820A: Cenko et al. 2006; GRB 060206: Wozniak 
et al. 2006; Monfardini et al. 2006), onset of the external 



shock (GRB 060206 and GRB 060210: Stanek et al. 2007; 
GRB 060418 and GRB 060607A: Molinari et al. 2007), and 
large angle emission (GRB 990123: Panaitescu & Kumar 2007). 
It is likely that the interplay of the many processes active after 
the GRB explosion may all affect the optical light curves, creat- 
ing the rich variety of observed behaviours. 

While GRBs with complex multi-peaked time profiles dis- 
playing no pulse width evolution with time seem to be explained 
more naturally through the mutual interaction of a wind of shells 
emitted at different times with different Lorentz factors (inter- 
nal shock model), a single FRED (fast rise exponential decay; 
Fishman & Meegan 1995) profile can still mat ch the expecta- 
tions of a single shell sweeping up the ISM (Fenimore et all 
fl996l) . In this scenario, the kinetic energy of a single ultrarel- 
ativistic shell is converted into internal energy of the ISM swept 
up; the shocked electrons radiate via synchrotron emission and 
inverse Compton scattering. In the simplest case of a thin shell 
ploughing into the ISM and emitting for a short time, the ex- 
pected time profile of the y-ray prompt emission is that of a sin- 
gle pulse with fast rise and slow decay. The cooling timescale 
of electrons is negligib le with respect to the hydrodynamical 
timescale in most cases dSari&PiranLll997l) . As a consequence, 
the rise time is determined by the emission time given by the 
crossing time of the shell by the reverse shock, while the decay 
time is dominated by the angular spreading timescale. 

In this paper we report on the y-ray, X-ray and optical ob- 
servations of GRB 070311, whose time profile is typical of a 
FRED. In particular, we focus on some properties shared by the 
y-ray and optical prompt emission and the late optical and X-ray 
afterglow in the light of the refreshed shock scenario. 

The paper is organised as follows: in Sect. [2] and Sect. [3] we 
describe the observations and the data reduction and analysis, re- 
spectively. Multiwavelength timing and spectral analysis of both 
the prompt and the afterglow emission is presented in Sect. [4] In 
Sect. [5]we discuss our results. Finally, in Sect. [6] we summarise 
our findings and conclusions. 

Throughout the paper, times are given relative to the on- 
set time of the GRB, which corresponds to 45 s prior to the 
INTEGRAL/IBAS trigger time, and the convention F(v, t) oc 

r a has been followed, where the energy index [} is related 
to the photon index F = /3 + 1 . 

All the quoted errors are given at 90% confidence level for 
one interesting parameter (Ay 2 =2.706), unless stated otherwise. 

2. Observations 

GRB 070311 triggered the INTEGRAL/IBAS in IBIS/ISGRI 
data on 2007 March 1 1 at 01 :52:50 UT and it was localised at RA 
= 5" 50 m 09:86, Dec = +0 3° 22' 29V3, with an error radius of 
2:5 dMereghetti et al.Ll2007l) . At the time of the burst Swift/BAT 
was pointing in almost the opposite direction and it would have 
been in the field of view two minutes later. The corresponding 
flux at the position of the burst in a 300-s image beginning 130 s 
after the INTEGRAL trigger is negligible. 

The y-ray prompt emission in the 20-200 keV energy band 
lasted about 50 s with a peak flux of 0.9 ph cirT 2 s _1 (1 s 
integration time) and a fluence of (2 7 3) x 10 -6 erg cm" 2 
dMereghetti et all 120071; ISazonov et aUl2007l) . 

The Rapid Eye Mounfl (REM; Zerbi et al. 2001) telescope 
reacted promptly and began observing 55 s after the GRB trigger 
time (see Sect. 13. It and discovered a bright fading IR counterpart 
within the INTEGRAL error circle at RA = 05 h 50 m 08 s .21, Dec. 
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= +03° 22' 3073 (J2000; Covino et al. 2007), corresponding to 
Galactic coordinates (/, b) of (202?7 66, — 1 1?998). The afterglow 
was soon confirmed by PAIRITEL dBloomll2007l) . 

The Swift narrow field instruments, XRT and UVOT, began 
observing at 7004 s after the trigger time. The XRT found an 
uncatalogued fading source at the position RA = 05 h 50 m 08*43, 
Dec. = +03° 22' 3070 (J2000), with an error radius of 378 and 
373 from the optical counterpart dGuidorzi et alll2007al) . No op- 
tical source was detected in correspondence of the optical and 
X-ray afterglows by UVOT down to V = 19.5 and B = 20.5 
(3cr) from 197-s exposures dHolland et all 120071) . 

The 2.2-m telescope of Calar Alto (CAHA) equipped with 
the Calar Alto Faint Object Spectrograph (CAFOS) observed the 
afterglow in R filter at 0.72 and 1.73 days after the burst. 

Because of the low Galactic latitude, the Galactic redden- 
ing along the directi on to the GRB is large: Eb-v = 0.763 mag 
(Schlegelet al., 1998). The Galactic extinction in each filter has 
been estimated through the NASA/IPAC Extragalactic Database 
extinction calculator. Specifically, the extinctio n in the other 
filters is derived through the parametrisation by ICardelli et alj 
(119891) : ;W = 2.53, A« = 2.04, A/ = 1.48, Aj = 0.69, Ah = 0.44, 
A K = 0.28 mag. 

3. Data reduction and analysis 

3.1. Gamma-ray data 

Figure [TJ shows the 18-200 keV background-subtracted 
time profile of GRB 070311 as recorded by INTEGRAL 
dSazonov et al ., 2007). The onset of the GRB appears to occur 
about 45 s before the trigger time. Hereafter, time will be re- 
ferred to the onset time, i.e. 45 s prior to the trigger time. From 
75 to ~ 125 s the signal drops below the sensitivity of the instru- 
ment. Interestingly, between ~ 125 and 225 s there is a hint of 
the presence of a faint gamma-ray tail. Although this feature is 
detected at ~ 2.5-cr confidence level (see thick cross in Fig. [TJ 
and so should be regarded as tentatively detected, it is interest- 
ing to see how it fits into the overall picture together with NIR, 
optical and X-ray observations of GRB 0703 1 1 . We address this 
issue below. 

The integrated spectrum in the 18-300 keV band is well fit 
by a power law with a photon index of Y y = 1.3+0.1 and 
a fluence of (3.0 ± 0.5) x 10~ 6 erg cm" 2 (20-200 keV). The 
spectrum shows no statistically significant high-energy cutoff. 
The peak energy lies above 80 keV during the prompt emission. 
There is also an indication of spectral softening in the course 
of the burst, with T y evolving from 0.8 + 0.2 during the rise 
phase to 1.45 ± 0-15 during the peak and decay of the emission 
dSazonov et aill2007l) . 

The conversion factor from rate to flux units in the 18- 
200 keV band is (1.2 + 0.2) x 10 -9 erg crrT 2 count -1 . 

3.2. Infrared/optical data 

Early time optical and near infrared (NIR) data were collected 
using the 60-cm robotic telescope REM located at the ESO 
La Silla observatory (Chile). The REM focal instruments con- 
sist of a NIR camera (REMIR), operating in the range 1.0- 
2.3 fim (z'JHK'), and an optical imager (REM Optical Slitless 
Spectrograph, ROSS) with spectroscopic (slitless) and photo- 
metric capabilities {VRI). A dichroic allows simultaneous ob- 
servations at optical and NIR wavelengths in two selected filters 

2 http://nedwww.ipac.caltech.edu/forms/calculator.html 
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Fig. 1. Histogram of the INTEGRAL 18-200 keV background- 
subtracted y-ray profile of GRB 07031 1 (integration time of 4 s; 
left axis). The thick cross shows a grouped bin of the y-ray pro- 
file, ~ 2.5 cr above the background. REM H (filled circles) and 
R (empty squares) magnitudes of the NIR/optical afterglow are 
also reported on the right axis. 



(for further information on REM and its capabilities, see Covino 
et al. 2004 and references therein). 

REM reacted promptly to the INTEGRAL GCN alert and 
began observing the field of GRB 0703 11 55 s after the burst 
trigger (36 s after the reception of the alert), following the event 
for ~ 1 hour. 

For the first ~ 500 s the REMIR observations have been per- 
formed using only the H filter with increasing exposure times, 
then all the NIR filters have been used in rotation. A similar ob- 
serving strategy has been adopted in the optical. R band observa- 
tions lasted ~ 1400 s for a total of 40 consecutive images. During 
the following ~ 2700 s, VRI images have been acquired in ro- 
tation but the optical transient was already below the instrument 
detection limits for the V and / filters. 

For both optical and NIR data sets, the reduction and the 
analysis followed standard procedures. The photometric calibra- 
tion for the NIR was accomplished by applying average mag- 
nitude shifts computed using the 2MASS3catalogue. The opti- 
cal data were calibrated using instrumental zero points, checked 
with observatio ns of standard stars in the SA95 Landolt field 
dLandold [19921) . or with the SDSS0 in the case of the zf filter. 

Figure [TJ shows the REM HR prompt measurements (filled 
circles and empty squares, respectively) together with the y- 
ray time profile, while Figure [2] shows the KHJz'R light curves 
(empty circles, crosses, filled squares, empty diamonds, empty 
squares, respectively). 

Further observations were acquired with CAFOS. This 
instrument is a focal reducer which allows direct imaging, 
spectroscopy and polarimetry. The detector used is a SITe 
2048x2048 pixel CCD providing a scale of 0753 pixel -1 and 
a circular field of view of 16' in diameter. In order to reduce the 
CCD readout time the observations of GRB 0703 1 1 were car- 
ried out trimming the CCD to a 1024x1024 window. The pho- 
tometric calibration was performed observing the standard field 
PG0942 dLandoltl 119921) at a similar airmass as the GRB field. 
Data reduction and analysis were carried out following s t andar d 
recipes by using ESO-Eclipse (v5.0) tools dDevillardl |200ll) . 

3 http://www.ipac.caltech.edu/2mass/ 

4 http://www.sdss.org 
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Fig. 2. Panchromatic light curves of the prompt as well as of the afterglow emission of GRB 070311: the flux in the 18-200 keV 
band by INTEGRAL (empty upside down triangles) and in the 0.3-10 keV band of the late afterglow by Swift/XRT (filled circles) 
are given on the right axis. Filled diamonds show the y-ray flux extrapolated to the 0.3-10 keV band. Magnitudes in J (filled 
squares), H (crosses), K (empty circles), z' (empty diamond) and R (empty squares) filters by REM and R by CAFOS (filled upside 
down triangles) are given on the left axis. A single H point from PAIRITEL (Bloom et al. 2007; empty triangle) and other R points 
from GCN circulars (Cenko 2007; Dai et al. 2007; Garnavich et al. 2007; Halpern & Armstrong 2007a,b,c,d; Jelmek & Prouza 
2007; Kann 2007; Wren et al. 2007; filled triangles) are shown as well. Upper limits are at 3-cr confidence level. The solid (dashed) 
line shows the best-fit power law (a = 1.06 + 0.08) and two pulses superposed on the H (early R, i.e. t < 10 4 s) filter curve. The 
dotted line shows the same power-law component renormalised to the first part (f < 10 5 s) of the X-ray curve (a x — an — ««)• The 
two earliest vertical dashed lines show the time interval used to extract an early SED (Fig. [Tol l, while the third vertical line shows 
the epoch of a late SED (Fig. [TTJi. See text. 



Aperture photometry was obtained with the GAIAEI (v3.2.0) 
package. REM and CAFOS photometry is reported in Table [TJ 

Hereafter, the magnitudes shown are not corrected for 
Galactic extinction, whilst fluxes as well as all the best-fit mod- 
els are. When the models are plotted together with magnitudes, 
the correction for Galactic extinction is removed from the mod- 
els. 

First we fitted the H light curve, which is the best monitored 
at early times. We tried to fit the points up to ~ 300 s with a sim- 
ple power law to test whether the fluctuations visible in Fig. [TJ 
might be just statistical around a power-law decay. The result- 
ing ^ 2 /dof is 17.5/7 with a chance probability of 1.4% and the 
power-law index turned out to be an = 0.5 ± 0.4 (1 cr). This 
probability gets even smaller if we combine the R and H points, 
as the latter trace the former although less significantly. We infer 
that the probability that the first 300 s NIR/optical profile is the 
result of statistical fluctuations around a simple power-law decay 
is lower than 1.4%. 

This motivated us to fit the early data by means of a more de- 
tailed model. The point from PAIRITEL at t = 1 160 s (iBlooml 
120071) was included in the H data set for this analysis. The fit 
shown in Fig. [2] (solid line) is the result of a power-law com- 
ponent with decay index of an = 1.06 + 0.08 with two FRED- 



5 http://star-www.dur.ac.uk/~pdraper/gaia/gaia.html 



shaped pulses overimposed, peaking at 119 and 180 s, respec- 
tively. The dashed line shows the power law fitting the early part 
of the R curve (f < 10 4 s) by fixing the index or = an and let- 
ting only the normalisation free to vary. Here we note the less 
pronounced enhancement in the R flux at the time of the pulses 
seen in H. A more detailed discussion of it as well as of the 
fitting models is reported in Sect. 14.11 

For both the K and J magnitudes there is one single detection 
and the remaining are upper limits. The first J upper limit at 
810 s preceding the detection at 1537 s is however inconsistent 
with the assumption aj — an extrapolated at the time of the 
two measurements. Comparing all the optical-NIR filters, a flux 
increase characterised by some variability seems to appear even 
after the first two optical peaks, i.e. after ~ 600 s. This, combined 
with the detections in J and z', could be a hint of a third flare. 

3.3. X-ray data 

The XRT began observing GRB 070311 on 2007 March 11 at 
03:49:34 UT, 7004 s after the INTEGRAL trigger and ended 
on 2007 March 19 at 22:38:54 UT, with a total net exposure 
of 81.8 ks in photon counting (PC) mode spread over 8.8 days. 
Table |2]reports the log of the XRT observations. 

The XRT data were processed using the FTOOLS software 
package (v. 6.1) distributed within HEASOFT. We ran the task 
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Table 1. REM KHJz'R and CAFOS R photometry of GRB 07031 1. 



Start Time 


End Time 


Exposure 


Mag a ' b 


Instr/Filter 


(s) 


(s) 


(s) 






554 


704 


150 


15.20 + 0.24 


REM/A" 


1281 


1431 


150 


>15.10 


REM/ A" 


59 


69 


10 


14.36 ±0.23 


REM/// 


75 


85 


10 


13.83 ±0.14 


REM/// 


90 


100 


10 


14.29 ± 0.20 


REM/// 


105 


115 


10 


14.21 ±0.19 


REM/// 


120 


130 


10 


14.70 ± 0.30 


REM/// 


138 


148 


10 


13.88 ±0.15 


REM/// 


154 


164 


10 


14.50 ± 0.24 


REM/// 


168 


178 


10 


14.37 ±0.22 


REM/// 


184 


210 


26 


14.58 ±0.19 


REM/// 


224 


274 


50 


15.38 ± 0.20 


REM/// 


303 


353 


50 


15.48 ±0.22 


REM/// 


372 


522 


150 


15.77 ±0.15 


REM/// 


1099 


1249 


150 


>16.4 


REM/// 


1826 


1976 


150 


>16.4 


REM/// 


2536 


2836 


300 


>17.3 


REM/// 


735 


885 


150 


>17.10 


REM/7 


1462 


1612 


150 


17.35 ±0.31 


REM/7 


2188 


2338 


150 


>17.10 


REM/7 


3198 


3498 


300 


> 17.24 


REM/7 


917 


1067 


150 


16.1 ±0.3 


REM/z' 


59.2 


89.2 


30 


17.49 ±0.22 


REM//? 


105.6 


135.6 


30 


17.86 ±0.21 


REM//? 


151.8 


181.8 


30 


17.73 ±0.19 


REM//? 


197.9 


227.9 


30 


18.12 ±0.27 


REM//? 


244.1 


274.1 


30 


18.19 ±0.35 


REM//? 


290.3 


320.3 


30 


18.31 ±0.35 


REM//? 


336.5 


597.0 


180 


18.91 ± 0.24 


REM//? 


613.0 


643.0 


30 


18.32 ±0.34 


REM//? 


659.2 


4295.6 


1530 


21.20 ±0.50 


REM//? 


62174 


74783 


8400 


22.57 ±0.12 


CAFOS//? 


149180 


160924 


10810 


22.60 ± 0.08 


CAFOS//? 



a Values are not corrected for Galactic extinction. 

b Errors at the 68% confidence level and upper limits (3 cr) are given. 



Table 2. XRT observation log of GRB 0703 1 1 . 



Sequence 


Obs Mode 


Start Time (UT) 


End Time (UT) 


Exposure 


Start Time a 


End Time" 










(s) 


(s) 


(s) 


00020052001 


PC 


2007-03-11 03:49:34 


2007-03-11 17:08:03 


18780 


7004 


54913 


00020052002 


PC 


2007-03-13 00:49:40 


2007-03-13 13:46:56 


9697 


169010 


215646 


00020052003 


PC 


2007-03-14 00:49:36 


2007-03-14 13:53:56 


9465 


255406 


302466 


00020052004 


PC 


2007-03-15 00:47:36 


2007-03-15 14:01:54 


9954 


341686 


389344 


00020052005 
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Since INTEGRAL trigger time. 



xrtpipeline (v.0.10.4) applying calibration and standard filtering 
and screening criteria. Data were acquired only in PC mode due 
to the faintness of the source. Events with grades 0-12 were se- 
lected. The XRT analysis has been performed in the 0.3-10 keV 
energy band. 



3.3.1. Temporal analysis 

Source photons were extracted from a circular region with a 20- 
pixel radius (1 pixel = 2'.'36; Fig. [3} and PSF-renormalised. The 
background was estimated from a four-circle region with a total 
area of 11.5 x 10 3 pixel 2 away from any source present in the 
field. When the count rate dropped below ~ 10~ 2 counts s _1 , we 
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Fig. 3. XRT image of the field of GRB 070311 obtained from 
82 ks PC data. The large circle shows the INTEGRAL position 
with an error radius of 2'.5 (90% CL). The cross shows the opti- 
cal afterglow position discovered by REM. The small circle is a 
20-pixel radius region centred on the XRT afterglow. 



made use of XIMAGE with the tool SOSTA, which corrects for vi- 
gnetting, exposure variations and PSF losses within an optimised 
box, using the same background region. 

The resulting 0.3-10 keV light curve is shown in Fig. [4] It 
was binned so as to achieve a minimum signal to noise ratio 
(SNR) of 3, except for the two last bins (with a SNR of 2.8 and 
2.6, respectively) as well as a minimum of 20 total counts. 

The light curve has been fit firstly with a double broken 
power law (solid line in Fig. @}, whose best fit parameters are 
the following: a x>x = 1.4 + 0.1, or x , 2 = -0.1^, Of x , 3 = 3.1+°-*, 
f x ,bi = 31±4ks, f Xjb2 = 1.7+° ' xlO 5 s Or 2 /dof = 16.8/24), where 
«x,i (i = 1,2,3) are the canonical initial, mid, and final decay 
slopes and f x b, (i = 1 , 2) are the two break times, respectively. 
The last point lies on the extrapolation of the initial power-law 
decay (dashed line in Fig.|4|i. From the X-ray data alone it is not 
possible to determine whether the shallow and final steep decay 
phases are the result of a late rebrightening after which the decay 
will resume to the pre-break behaviour. 

We also adopted two alternative models: the combination 
of a FRED-shaped pulse with a power law (dash-dotted line 
in Fig. |4j a more detailed descript i on fol lows in Sect. 14.1.21 ), 
and the model by Willingal e et al.l (12007). This model is the 
combination of two components, the prompt and the afterglow 
(according to the terminology introduced by these authors), de- 
scribed with the same functional form, which consists of a com- 
bination of an exponentia l and a power law. Following the no- 
tation bv lWUHngaleetaTI d2007h . we fixed the rise time of the 
prompt component to the onset time: f p = 0. We also fixed the 
time when the power-law component of the prompt takes over, 



T p = 5 x 10 3 s, i.e. prior to the beginning of the X-ray observa- 
tions, when the X-ray decay is already dominated by the power 
law. The best-fit model is shown in Fig. [5] The best-fit parame- 
ters are the following: a p = 1.6 ± 0.2, f a = (1.0 + 0.4) x 10 5 s, 
T R = (2.9 + 1.0) x 10 5 s, a a = 3.5 + 0.7 ( X 2 /dof = 24.9/24). 
We point out that within this model the late rebrightening corre- 
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Fig. 4. Top panel: X-ray afterglow light curve of GRB 070311 
obtained with XRT in the 0.3-10 keV energy band. The solid 
line shows the best fit obtained with a double broken power 
law, while the dashed line shows the extrapolation of the ini- 
tial power-law decay. The dashed-dotted line shows the best- 
fit combination of a FRED-like pulse superposed on a power 
law. Labelled vertical slices correspond to four different regions 
where spectra have been extracted. Mid panel: fractional residu- 
als with respect to the double broken power-law model. Bottom 
panel: photon index (error bars are 1 cr). 




Time since GRB onset (=trigger time - 45 s) [s] 

Fig. 5. Top panel: X-ray afterglow light curve of GRB 0703 1 1 
obtained with XRT in the 0.3-10 keV energy band. The curve 
has be en fit with the two-component model by Willingale et al.l 
(2007): the prompt (dashed), the afterglow (dotted) and their 
combination (solid). Labelled intervals are the same as in Fig. [4] 
Bottom panel: fractional residuals with respect to the model. 



sponds to the onset of the afterglow component: this forced us 
to decouple f a from T p and trea t the f ormer as a free parame- 
ter, unlike what I Willingale et al.l (|2007|) did for all of the GRBs 
of their sample. We las t note that the two-component model by 
Willingal e et al.l (120071) usually accounts for both the prompt 
and the afterglow emission. In our fit, we have only consid- 
ered the late-time light curve. The X-ray rebrightening observed 
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Fig. 6. Top Panel: X-ray afterglow photon spectra of 
GRB 070311 obtained with XRT in the 0.3-10 keV en- 
ergy band corresponding to the temporal interval A (crosses) 
and BCD (circles), respectively. Solid lines show the corre- 
sponding best-fit absorbed power laws. Bottom Panel: residuals 
with respect to the corresponding best-fit models. 



as high as that derived from the radio, interpreted as due to 
the presence of molecular gas ( Arabadiis & Bregman, 1999; 
iBaumgartner & Mushotzkvl [2006). Therefore we conclude that 
the A^h we measure from the X-ray spectrum is consistent with 
the Galactic value expected in the direction of GRB 07031 1, al- 
beit we cannot exclude some intrinsic absorption. 

Comparing the best-fit parameters obtained for the different 
spectra, while the absorption seems constant, we find the sugges- 
tion (2.2(f) for a softening of the photon index F x , from 2.05*2^ 
(interval A) to 2.5 + 0.2 (intervals B, C, D merged together). 
The corresponding spectra with the best-fit models are shown in 
Fig. [6] 

We note that the change of the spectral index, A/3 X = Ar x , is 
consistent with the canonical value of 1 /2 expected in the stan- 
dard synchrotron model when the cooling frequency crosses the 
observed passband (X rays, in this case). In this case, this pas- 
sage would have occurred approximately between intervals A 
and BCD (Fig. H. 

The excess visible between 4 and 5 keV in the residuals of 
the spectrum with respect to the absorbed power law appears to 
be ~ 2.5-cr significant after properly rebinning and in the most 
favourable case (Pig. 



in GRB 0703 1 1 would therefore constitute a third component, 
which is only rarely seen in afterglow light curves. 

3.3.2. Spectral analysis 

In order to detect possible spectral variations connected with 
changes in the X-ray light curve, we extracted the 0.3-10 keV 
spectrum in four different time intervals, labelled "A", "B", "C" 
and "D", corresponding to the initial steep decay, the beginning 
of the shallow phase, the end of it (or the peak of the rebright- 
ening, according to the alternative description of the X-ray light 
curve discussed in Sect. 13. 3. Ti l, and the final steep decay, respec- 
tively (see Fig.|4|. 

Source and background spectra were extracted from the 
same regions as the ones used for the light curve (Sect. 13.3.11 ), 
except for intervals "B" and "C" separately. For these intervals, 
due to their poor statistics, it was not possible to group the en- 
ergy channels so as to have a Gaussian distribution of the num- 
ber of photons per grouped channel. H ence, f or bot h spectra we 
replaced the x 1 with the C statistics (Cash|, 119791) . which has 
proven to be useful whenever the Gaussian approximation does 
not hold (e.g. when the number of photons per bin is less than 
10), provided that the contamination of background photons is 
negligible. In order to ensure this, for both spectra "B" and "C", 
a 10-pixel radius circular region was used. 

The ancillary response files were generated using the task 
xrtrakarf . Spectral channels were grouped so as to have at least 
20 counts per bin, except for "B" and "C". Spectral fitting was 
performed with xspec (v. 1 1.3.2). 

All the spectra can be modelled with an absorbed power law 
with the xspec models wabs pow. We assumed th e photoelectric 
cross section by lMorrison & McCammonl (11983). Results of the 
best-fit parameters are reported in Table [3] 

The Galactic neutral Hydrogen column density along the 



GRB direction from 21 -cm-line radio surveys is N, 



(Gal) 
HI 



2.3 x 



10 21 cirT 2 dKalberla et all 120051) . The X-ray absorption found 



from spectral fitting is significantly higher than A^ dl) : specif- 
ically, it is about twice as high. We do not interpret this as 
evidence for intrinsic absorption: for low Galactic latitudes 
(b < 25°) the A^h measured from X rays is about twice 



4. Multi-band combined analysis 

4.1. Light curves fitting 

We analysed the different decays in different energy bands as 
homogeneously as possible. Motivated by the possible interpre- 
tation of the FRED in the external shock context, and by the 
shape of the //-filter pulses similar to that of the prompt emis- 
sion in the 18-200 keV band, which looks like a typical FRED, 
we first fitted the 1 8-200 keV pulse adopting the profile used by 
iNorris et al.ldl996l) : 



F(t) 



A exp 
A exp 



(^max t V 
<T r / 

(* ^max \ 
cr d / 



, t < fmax 
, t > ? m ax 



(i) 



where f max is the peak time, cr T and o-& are the rise and decay 
times, respectively, A is the normalisation and v is the peaked- 
ness (when v = 1 the profile is a simple exponential, when v - 2 
it is a Gaussian). The best-fit parameters for the y-ray pulse are: 
fpeak = 39.0 + 0.8 s, a- t = 8.5 ± 1.0 s, <x d = 15.0 +1.7 s. The 
peakedness was found to be v = 1.08+0.13 starting with an ini- 
tial value of v = 1 . The result is shown by the solid line in Fig. [7] 
We calculated the flux density at 88 keV, which corresponds to 
the energy at which the flux density equals the average flux den- 
sity in the 18-200 keV band assuming B y = 0.3. 

To fit the curves, we adopted two alternative approaches. 
First we focused on the early pulses and fitted them in terms of 
two overlapping FRED-shaped pulses (pulses model). We chose 
to model the shape of the pulses with a FRED, first because 
this fits well, and, secondly, because it allows a better compar- 
ison with the results of the fit of the y-ray pulse. Alternatively, 
we added a power-law continuum, in the assumption that the 
afterglow contribution is not negligible at this time (pl+pulses 
model). Both H and R profiles have been corrected for Galactic 
extinction. Given the less dense sampling of the H and R curve 
with respect to the y-ray one, the peakedness was fixed to the 
best-fit value reported above. 

Finally we fitted the late rebrightening seen in R and X-ray 
with a single pulse superposed to a power law (pl+pulse model). 
Best-fit parameters of the models are reported in Table [4] 
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Table 3. Best-fit parameters of the 0.3-10 keV spectrum of the X-ray afterglow. The model is an absorbed power law (xspec model: 
wabs pow). 



Interval Start time 


Stop time Nh 


r x 


Mean flux 


Af 2 Mof 


(s) 


(s) (10 21 cm" 2 ) 




(10~ 13 erg cnr 2 s _1 ) 





11.0/16 
351.5 (63.4%)' 
320.4 (58.2%)' 
12.4/9 
10.3/13 
6.2/7 



" Cash statistics (C-stat; Cash 1979) and percentage of Monte Carlo realisations that had statistic < C-stat. We performed 10 4 simulations. 
Photons were extracted from a circular region with a 10-pixel radius. 
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765966 
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5.5! 2 | 2.5 + 0.2 1.9 + 0.1 

4.4+j' 2.4 + 0.4 1.1+0.4 
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Time since GRB onset (=trigger time - 45 s) [s] 

Fig. 7. INTEGRAL background-subtracted y-ray flux density at 
88 keV (crosses; integration time of 4 s) magnified 30 times, 
REM H (filled circles) and R (empty squares; magnified 4 
times) flux densities of the NIR/optical afterglow (corrected for 
Galactic extinction). Solid, dashed and dotted lines show the 
best-fit models of the y-ray, H- and /?-band pulses, respectively. 
The thick cross is the same as that in Fig.[T] 



4.1 .1 . pulses model 




100 120 140 160 180 200 220 240 260 280 
Time since GRB onset (=trigger time - 45 s) [s] 

Fig. 8. Early R optical flux (dashed squares) compared with the 
time-integrated flux expected from the best-fit model of the H 
curve rescaled by a factor of 0.55 + 0.06, which minimises the 

x 1 - 



fitting the H pulses is consistent with the simultaneous R mea- 
surements. 

R points are displayed in Fig. [7] with empty squares and the 
best-fit model is represented with the dotted line, all magnified 
by a factor of 4 for the sake of clarity. 



The first H pulse peaks 80 s after the y-ray peak. The pulse shape 
is different: the rise time of the H pulse is about twice as long, 
while its decay time is about 4 times longer. The flux density 
normalisation constant, Ahi = 4.9 + 0.3 mJy, is about 17 times 
that of the y rays at 88 keV. The second H pulse follows the first 
by ~ 62 s. Because of the poor sampling of the rise, only the 
decay is better constrained and turns out to be about twice as 
long as the decay time of the first pulse, while the normalisation 
of the second pulse is about half as big: Am - Am/2, so that the 
fluence during the decay is similar to that of the first pulse. The 
dashed line in Fig.|7]shows the sum of both pulses fitting the H 
points (filled circles). 

To evaluate whether the R points simultaneous to the promi- 
nent pulses seen in H are consistent with being derived from the 
same profile as H within uncertainties, we integrated the best- 
fit model of the H curve over the R time bins. We determined 
the rescaling factor / Des t that minimises the x 1 between the mea- 
sured R and the expected rescaled mean H flux. The result is 
shown in Fig. [8] and corresponds to /best = 0.55 + 0.06, yield- 
ing ^ 2 /dof = 0.36/3. We conclude that the time profile best- 



4.1.2. pl+pulses model 

Alternatively to the description of Sect. I4.1.1I here we assume 
that the contribution of the power-law continuum is not negligi- 
ble since the very beginning of the NIR/optical observations. The 
result is shown in Fig. [2] For the same reasons as in Sect. I4.1.T1 
we first fitted the H profile. The power-law continuum turned 
out to have a slope of an = 1.06 + 0.08 (solid line in Fig. [2j. 
Then we added the same combination of pulses as that found 
in Sect. 14.1.11 by allowing the single normalisations and releas- 
ing gradually some of the parameters. Table|4]reports the best-fit 
values. Similarly to Sect. 14.1.11 the fit for the R filter was done 
by allowing just a scaling factor of the H profile. This turned out 
to be /best = 0.55 + 0.07, i.e. the same as that obtained with the 
previous model (Sect. 14.1.11 ). The R best-fit profile is shown in 
Fig.|2]with a dashed line. The^ 2 /dof of the best-fit model are ac- 
ceptable: 6.0/7 and 12.1/8 for the H and R profiles, respectively. 

The main differences from the results obtained in Sect. |4.1.T| 
concern the shorter rise times in this case, that were fixed to 
0.1 s. This is a consequence of having increased the continuum 
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Fig. 9. Late time rebrightening and correspondent best-fit mod- 
els (power law and a pulse) seen in R filter (triangles, thick 
dashed line) and 0.3-10 keV band (circles, thick solid line). 
Thin lines show the correspondent power-law component alone. 
X-rays: the dotted line shows the rescaled version of the y-ray 
prompt pulse combined with the underlying power law, while 
the dashed-dotted line represents the sum of two power laws. 



component here represented by the power law. The pulse shape 
is less constrained, for the same reason. The second pulse still 
has a longer decay than the first, while its peak intensity is now 
comparable with that of the first. 

The dotted line in Fig. [2] shows the best-fit power law of the 
X-ray data up to 10 5 s obtained by fixing a x — an — o:r. The 
fit is acceptable: ;^ 2 /dof = 24.3/21 (chance probability of 28%). 
This description of the X-ray light curve is alternative to that 
given in Sect. 13.3.11 

We tentatively extrapolated the y-ray flux to the 0.3-10 keV 
band, assuming the simple power law with T y = 1.3 + 0.1 
(Sect. I3.lt . This assumption relies upon the fact that the peak 
energy is likely to lie above 80 keV because of the hardness of 
the photon index T y . We also applied the suppression factor due 
to the soft X-ray absorption measured in the X-ray spectrum. 
The result is shown with filled diamonds in Fig. [2] Interestingly 
the single power law derived from fitting the early H and R con- 
tinuum components, a x = 1 .06, matches the level of the extrap- 
olated prompt emission flux in the X-ray band. 

4.1.3. Late tf/X-ray rebrightening 

From Fig. [2] we note that the power-law continuum, over which 
the late (f ~ 2 x 10 5 s) R rebrightening sets in, looks brighter 
than expected from extrapolating the early best-fit model by ~ 1 
order of magnitude. For this reason we fitted it separately from 
the early R light curve. We adopted a combination of a power- 
law and a single pulse (Eq. [1). The power-law index was fixed 
to the value found in Sect. 14.1.21 for the early part of the R, i.e. 
an = 1.06. The result is shown in Fig.|9](thick dashed line). We 
did the same for the contemporaneous 0.3-10 keV profile (thick 
solid line in Fig. [9j. We initially fixed a x = or in analogy with 
Sect. 14.1.21 but we had to release it because of the badness of 
the fit and it turned out to be a x = 1.47 + 0.20, more consistent 
with the broken-power-law fit given in Sect. 13.3.11 The X-ray 
coverage of the late rebrightening is not as detailed in catching 
the peak time as for the R filter, nonetheless from Fig.|9]we can 
confidently affirm that X-rays peak earlier than optical, which is 
rising or flattening during the beginning of the monitored X-ray 
decay. Best-fit parameters are reported at the bottom of Table |4] 



The dashed-dotted line in Fig. [9] represents the case when 
the X-ray late rebrightening is fit with the sum of two power 
laws, the second of which turns on between 6 x 10 4 and 2 X 
10 5 s and models the rebrightening superposed to the first power 



law. The best-fit indices are a 



(pp) 



1.4 + 0.1 and a. 



(pp) 



3.5 ± 



x,l * ■ x,2 

0.5, respectively (1 cr confidence). For both components the time 
origin was fixed to the GRB onset time. The time at which the 
second power law sets in cannot be estimated from our data; 
however, this is irrelevant for determining the power-law indices. 

Compared to the R band, the poorer coverage of the X-ray 
peak reflects upon bigger uncertainties on the best-fit parame- 
ters of the X-ray profile. Interestingly, rise and decay times are 
similar and the ratios o-<i/o- r are 2.3 and 2.6 for the R filter and 
X-ray, respectively, i.e . the pulses resem ble the typical shape of 
a prompt y-ray pulse (Norri sTt al.L 1 19961) . The peak intensities 
are 6.5 ± 0.9 and 18+9 times the value of the correspondent un- 
derlying power-law at peak time for the R filter and X-ray band, 
respectively. The X-my/R time lag amounts to (9.1 +3.3)x 10 4 s. 
More simply, the late R rebrightening peaks twice as late as the 
X-ray profile. A word of caution is needed about the evaluation 
of this temporal lag: the uncertainty could be larger than our es- 
timate, which is constrained by the choice of the functional form 
used for fitting. Nevertheless, the evidence for a positive lag is 
apparent, regardless of the fits. 

4.1 .4. Prompt y pulse/late X-ray rebrightening 

Notably, the ratios between the temporal parameters best fitting 
the pulse of the y-ray prompt profile (Sect. I4.1.U and their cor- 
respondent ones fitting the late X-ray rebrightening (Sect. I4.1.3t 
are all comparable: f pe ak,iatex/fpeak,y = 2600 + 800, cr r i ateX /cr^ = 
3400 ± 3000, cr djlateX /cr d y = 4400 ± 800. We tried to fit the late 
X-ray rebrightening with the combination of the same power law 
as that obtained in Sect. 14.1.31 and a rescaled version of the y-ray 

prompt pulse: f pea k,lateX = /s,X (f p eak,y - to), Cr r;lateX = /s,X CTt,y, 

c"d,iatex = /s,x o~&,y We left three parameters free to vary: the 
scaling factor / s x , the time origin to and the normalisation of 
the pulse. The choice of letting the time origin vary was moti- 
vated by the peak time of the y-ray pulse being very sensitive 
due to its smallness. The result is shown in Fig. [9] (dotted line). 
The best-fit parameters are the following: / s X = 5700 + 700 and 
t = -19.0 ± 2.4 s (x 2 /dof = 19.3/27). The potentially strong 
implications on the interpretation of this result are addressed in 
SectS 



4.2. SED evolution 

We derived two SEDs at different epochs. The early one com- 
prises the NIR pulses seen at the beginning of the REM follow- 
up observations and lasts from 104 to 273 s. This SED consists 
of detections in two filters, H and R, and a 3-cr upper limit of 
6.3 yuJy in the 1 8-200 keV band (see Fig. [2j». Given the high vari- 
ability of the H curve, especially when compared with that of R, 
the mean H flux was derived by integrating the best-fit profile of 
the H curve described in Sect. 14. l.TI over this time interval. The 
simultaneous mean R flux was calculated assuming the best-fit 
rescaling factor (/best = 0.55 + 0.06) between H and R derived in 
Sect. |4.1| for the same time interval. The fit with a simple power 
law yields /0NiR-o P t = 0.65 ± 0.21. The extrapolation of the fit to 
the y-ray band is consistent with the upper limit (Fig. [Tot. 

Interestingly, the spectral index is consistent with that at the 
end of the y-ray pulse (Sect. 13.11 ), fi 7 = 0.45 + 0.15, and might 
be suggestive of an unbroken power law segment from NIR to 
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Table 4. Best-fit parameters of the multi-band light curves of the afterglow of GRB 0703 1 1 . Uncertainties are 1 <x. Values of frozen 
parameters are reported among square brackets. 



Energy 
band/filter 


Model 


Component 


a 


t 

peak 

(s) 


(s) 


0"d 
(s) 


A 
(mjy) 


V 


18-200 keV 


PULSES 


pulse 


- 


39.0 ±0.8 


8.5 ± 1.0 


15.0 ± 1.7 


0.29 ± 0.05 b 


1.08 ±0.13 


H 


PULSES 


1 st pulse 


- 


119.0 ±2.2 


18.2 ±4.8 


63.6 ±5.2 


4.9 ±0.3 


[1.08] 


H 


PULSES 


2 nd pulse 


- 


180.5 ±9.4 


0.5 ± 3.6 


110 ± 15 


2.6 ± 0.2 


[1.08] 




PITT SFS 


1 s * mil 




n 1 9 oi 


H8 21 


[63.6] 


2.7 ± 0.2 


n 08i 


R 


PULSES 


2 nd pulse 


— 


[180.5] 


[0.5] 


[110] 


1.4 ±0.1 


[1.08] 


H 


PL+PULSES 


Pi 


1.06 ±0.08 


- 


- 


- 


3.6 ± 2.0 C 


- 


H 


PL+PULSES 


1 st pulse 


— 


119.0 ±0.1 


[0.1] 


27 ± 18 


2.2 ± 1.0 


[1.08] 




PT 4-PII1 SFS 


2 nd nulse 




181 1+06 


ro n 

LV. i J 


41 ± 17 


2.6 ± 0.8 


ri 08i 


R 


PL+PULSES 


Pi 


[1.06] 








2.0 ± l.l c 




R 


PL+PULSES 


1 st pulse 




[119.0] 


[0.1] 


[27] 


1.2 ±0.6 


[1.08] 


R 


PL+PULSES 


2 nd pulse 




[181.1] 


[0.1] 


[41] 


1.4 ±0.4 


[1.08] 


late/? 


PL+PULSE 


Pi 


[1.06] 








11.0±0.9 dc 




late R 
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2.85 ± 0.54 f 


7.4 ± 0.7 f 


36±4 C 


0.91 ±0.34 


0.3-10 keV 


PL+PULSE 


Pi 


1.47 ±0.20 








1.3 ±0.1 d g 




0.3-10 keV 


PL+PULSE 


pulse 




10±3 f 


2.9 ± 2.5 f 


6.6 ± 0.9 f 


24 ± 12 g 


0.89 ±0.56 



a Time since GRB onset (corresponding to 45 s before the trigger time). 

b Flux density corresponding at 88 keV. 

c Flux at t = 100 s. 

d Flux at t = 10 5 s. 

c Units of yujy. 

f Units of 10 4 s. 

g Units of 10~' 3 erg cm -2 s -1 . 



y at this epoch. However, the large value of Eb-v makes the 
uncertainty 5Ay on the Galactic extinction quite large, 5Ay > 
0.5 mag. Should the Galactic extinction in V be larger (smaller) 
by 0.5, the corrected spectral index of the NIR/optical data alone 
would bey8 N r R _ p t = 0.3+0.2 (^ N iR_ opt = 1.0+0.2). Thus, nothing 
conclusive can be said about the possible presence of extinction 
in excess of the Galactic one. 

The late SED (Fig. [TTT > was extracted at 2.6 x 10 5 s, cor- 
responding to the beginning of the decay following the late re- 
brightening in R simultaneously with the final steep X-ray decay 
(se e Fig. [2]). This SED includes a single R measurement taken 
by lHalpern & Armstrong! d2007dl) and the XRT spectrum BCD 
(Sect. 13.3.21 Table [3]) renormalised through the best-fit power- 
law segment of the X-ray light curve (Sect. 13.3. It at the epoch 
of the R point. The SED was fit with an SMC-extinguished (Pei 
1992; in the observer frame), X-ray photoelectrically absorbed 
broken power law with f3 x - jSnir-opi = 0.5. The A^h was fixed 
to the value already derived from the corresponding X-ray spec- 
trum, i.e. 5.5 x 10 21 crrT 2 , consistent with being entirely Galactic 
(Sect. |3.3.2| >. We found j8nir-opi = 1.0 ± 0.2, consistent with 
j6 x = 1.5 from the X-ray spectrum best fit (Table [3]). The break 
frequency turned out to be Vb = (3.0 + 0.9) x 10 17 Hz and 
A v = 0.80 ± 0.15 mag Of 2 /dof = 9.6/12). Figure \U\ shows the 
SED: the dashed (dotted) line represents the best-fit model with 
optical extinction and X-ray absorption taken out (shown). 

Alternatively, we tried to fit the SED with a single absorbed 
power law. If we left the slope free to vary, the fit is driven by 
the more numerous X-ray points, leading to f3 ~ 1.5 and Ay = 
6 ± 3 mag. However, if we correct the early SED for such an 
extinction, the intrinsic optical spectrum would be nonphysically 
blue (/i op t = -3). Otherwise if we impose a single power law 
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10 14 10 15 10 16 10 17 10 18 10 19 10 so 
Frequency (Hz) 

Fig. 10. NIR-optical/y-ray SED between 104 and 273 s after the 
GRB onset time. The dotted line shows the best-fit power law, 
with /0NiR-opt = 0.65 + 0.21. The upper limit is at 3 <x. 



comparable with that found in the case of a broken power law. 
However, although the fit cannot be rejected (y 2 /dof = 18/13), 
the residuals of the X-ray points with respect to the model show 
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10 14 10 15 10 16 10 17 10 18 10 19 
Frequency (Hz) 

Fig. 11. Optical/X-ray SED at 2.6 x 10 5 s, around the peak of 
the late rebrightening. The dashed (dotted) line shows the best- 
fit absorbed broken power law, corrected (not corrected) for the 
optical extinction and X-ray absorption. The Nh was fixed to 
5.5 x 10 21 cuT 2 found from fitting the X-ray spectrum alone, 
while Ay was found to be 0.80 + 0.15 mag, adopting an SMC 
profile (at z = 0). 



a trend. Therefore, a simple absorbed power law does not seem 
to be a good representation of our data. 

Because of the unknown redshift z of GRB 0703 1 1 , we cau- 
tion that the values of Ay computed in the fits above are calcu- 
lated for z — assuming a specific extinction law (SMC model), 
so they must be taken as indicative upper limits to the corre- 
sponding rest-frame values. Unlike for the early SED, the X-ray 
absorption in terms of Nh is known to be consistent with the 
expected Galactic value. Due to the unknown redshift, it is not 
possible to set an upper limit to the possible intrinsic rest-frame 
Nu, z absorption. For this reason and because of the upper limits 
on Ay, nothing can be inferred on the amount of dust and gas 
along the line of sight to the GRB progenitor. 



5. Discussion 

Together with the early optical/NIR flares possibly accompanied 
by a y-ray tail, the late, bright and long rebrightening seen in 
X-ray and R band probably is the most intriguing feature of 
GRB 0703 1 1 . Although such late brightenings or flares are not 
unprecedented , only a few bursts have exhibited them so far: e.g., 
GRB 970508 dPiro et all [19981: iGalama etjD, [1998J) the short 
GRB 050724 dCampana et all 1200 6; Males ani et all 120071). the 
z = 6.3 GRB 050904 dCusumano et all 120071: IWatson et ail 
2006). We also note the presence of the shallow decay phase 
preceding the X-ray brightening (or corresponding to its gradual 
onset). Due to the lack of contemporaneous coverage in the R 
band, we did not observe a similar shallow decay phase in the 
optical bands. However, from Fig. [2] we infer that it must have 
taken place (e.g., some energy injection between ~ 10 4 s and 
~ 10 5 s) to power the continuum component of the 7?-filter de- 



cay at late times, which is significantly above the extrapolation 
of the early data. 

In the case of GRB 0703 1 1 , we note that At/t ~ 1 for the late 
rebrightening in both X-ray and R filter, so timescales arguments 
cannot be used against an external shock origin for it. We note 
that this is in contrast to what has been observed for other GRBs 
(e.g. GRB 050502B; Burrows et al. 2005, Falcone et al. 2006; 
GRB 050724; Barthelmy et al. 2005, Campana et al. 2006), in 
which the rise and decay of the flares are too steep and require a 
resetting of the time origin. 

In addition, the FRED-like shape of the 18-200 keV prompt 
light curve and consequent lack of high variability, classi- 
cally used to argue against an external origin of the prompt 
emission of complex bursts (Sari & Piran 1997; but see also 
Dermer & Mitman 1999), indicates that the prompt emission of 
GRB 07031 1 is potentially consistent with it (e.g., see Kumar & 
Panaitescu 2003), although unlik ely due to synchrotron radiation 
dRamirez-Ruiz & Granotj , l2006l) . 

Motivated by this and by the analogies found between 
prompt and late afterglow, hereafter we try to interpret the ob- 
servations of GRB 0703 1 1 in the light of an external shock ori- 
gin. While the y-ray pulse is thought to be produced during the 
deceleration of the shell against the surrounding medium, we 
consider two alternative explanations for the late rebrightening: 
either a refreshed shock by a second shell emitted after the first 
one, or a single shell whose forward shock encounters a density 
bump. 

The description of the R and X-ray late rebrightening as 
a simple power law with a slope changing after the peak 
(Sect. |4.1.3I >, is consistent with the interpretation that the late 
rebrightening is due to a thin shell that caught up with the shock 
front of the blastwave at l ater times, as expected in the so-called 
refreshed-shock scenario (Rees & Meszaros, 1998). In this pic- 
ture, the fireball rebrightened and soon afterwards (~ 2 x 10 5 s) 
turned off: the decay is the result of two components: the pre- 
existing power law and the high-latitude radiation left over by 
the refreshed shock. This is the case when the energy of the im- 
pacting shell is lower than that of the fireball (as also suggested 
by the ratio between the radiated energy during the late flare and 
the prompt emission), so that no noticeable step in the power-law 
decay is observed, but still luminous enough to produce a peak 
on a short timescale. 

According to t he so-called "cu r vature ' ' effect 

( Fenimore et all 119961 ; iKumar & Panaitescul 120001; iDermerL 
|2004|) . the steepness of the fast decay phase of several X-ray 
afterglows as well as of the y-ray tail of single pulses is due to 
high-latitude radiation Doppler-boosted in the observer fixed 
energy band. A relation is expected between temporal and 
spectral index: a = /3 + 2 in its simplest form in the case of a 
thin shell emitting for a short time, where the time origin must 
correspond to the time of the GRB explosion as measured in the 
detector frame. 

Notably, for GRB 07031 1 both the decay of the y-ray pulse 
and the longer decays in X-rays and /^-filter of the late re- 
brightening are consistent with the high-latitude closure rela- 
tion and with the same time origin, i.e. the GRB onset time: 
a y ~ Py+2 — 2.3+0.1 and a x — fi x +2 = 3.5+0.2. This is proven 

by the measured slope of the second power law, aff = 3.5+0.5, 
which describes the flare decay, in addition to the underlying 



power law with a. 



(PP) 

x.l 



1.4 ± 0.1 (Sect. I4T31 . Therefore, if 



we adopt the combination of a power law and either a pulse or 
another power law for the late rebrightening (Sect. 14.1.3b . the 
slope of the pulse decay is still consistent with a curvature ef- 



12 



C. Guidorzi et al.: Prompt and afterglow study of GRB 07031 1 



feet: ar*f = fi x + 2 In addition, from Sect. 15X51 we know that 
if we move the reference time backward to 19.0 + 2.4 s before 
the GRB onset, the late brightening is well fit by a time-rescaled 
version of the prompt y-ray pulse. What is remarkable is that the 
peak time scales exactly in the same way as the pulse rise and 
decay times: Af yi p rompt /f r 

, prompt — Af x ,iate/f x ,iate- This suggests that 
either two shells generated the prompt pulse as well as the late 
rebrightening or, alternatively, that a unique shell caused both by 
encountering two density enhancements. We note that the late X- 
ray hump occurs later (few xlO 4 s) than what is more commonly 
seen in the typical flat phase (~ 10 3 s), also usually interpreted 
in terms of energy injection from the central engine. The detec- 
tion in the R filter rules out a very high redshift, so such a late 
refreshing may be either due to a more distant radius of emission 
of the forward shock (e.g. due to low density of the surrounding 
medium) or due to a slow shell catching up with the shock front. 

The 1 8-200 ke V fluence emitted in the prompt amounts to 
(3.0 + 0.5) x 10~ 6 erg cm -2 by integrating the best- fit model 
of Sect. 14.1.11 while the net 0.3-10 keV fluence of the late re- 
brightening is (2.4 + 1.2) x 10~ 7 erg cirT 2 , i.e. lower by one or- 
der of magnitude. If we extrapolate the prompt spectrum to the 
XRT band as in Sect. 14.1.21 and correct for the X-ray absorption, 
the fluence of the late X-ray rebrightening becomes comparable 
with that of the prompt extrapolated to the same energy band. 
Differently, the late R hump has a time-integrated flux about one 
order of magnitude larger than the early pulses seen in the same 
filter. The larger energy content of the late shock, with respect 
to the early one seen in optical, might be explained with an in- 
creasing efficiency in converting the kinetic energy of the blast- 
wave into radiation. This implies a change of the microphysical 
parameters of the afterglow, as suggested to explain the shallow 
decay phase of other GRBs dloka et all[2006tlPanaitescdl2006h . 
This can also be explained more simply as due to the fact that the 
vF v broadband spectrum peaks at lower energies at later times, 
so that the energy release in the observed R filter during the late 
rebrightening i s larger than durin g the prompt or soon after that. 

Following Iok a et alj d2005l) . we can find which scenarios 
may be compatible with the observed flux increase, AF/F * 10 
(both in R and X; Sect. 14.1.31 ), and the observed At/t ~ 1: 
from their fig. 1, the late hump of GRB 070311 is compatible 
with the refreshed sho ck scenario. Also the patchy shell model 
(Meszaros et al.l[l998l) . characterised by an anisotropic emitting 
surface of the fireball, is not ruled out, while the scenario of a 
density bump causing the late rebrightening is ruled out, un- 
less many clumps of matter are illuminated at the same time. 
Therefore, the bright flux of the rebrightening, AF/F as 10, 
seems to disfavour the interpretation of a density medium en- 
hancement, causing the observed flux increase. However we note 
that this is debated: e.g., according to Dermer & Mitman (1999) 
and Dermer (2007a), this could be produced by parts of the blast- 
wave in which most of the energy has not been converted into 
radiation, yet, while impacting on a dense and thick clump of 
matter. In this scenario, the same shell would be responsible for 
both the prompt and the late rebrightening: e.g., the high-latitude 
emission observed during the decay of the late hump might be 
the result of the blastwave finally reaching the wind-termination 
shock. However, the interpretation of the late rebrightening as 
due to a density bump seems disfavoured by the crossing of the 
X-ray band by the cooling break v c , whereas the observed fre- 
quency must be below v c . In addition to that, the remarkable 
flux enhancement observed, AF X /F X = 18 ± 9 (Sect. 14.1.31 ), 
makes this interpretation less probable. In both scenarios (re- 
freshed shock produced by another shell or density bump), the 
scaling factor of the timescale of the late brightening with re- 



spect to the early pulse, / s = 5700 + 700, could result from 
the interplay of two factors: the increase of the visible portion 
of the emitting surface and the fact that the blastwave Lorentz 
factor has decreased by a factor of Ti/r2 head on, thus stretch- 
ing the timescale by (H/I^) 2 . Therefore, since Ri < R2, from 
/ s = R0JR1 (H/rz) 2 , it must be H/r, < 4% = 75 ± 5. 

The early NIR/optical pulses detected soon after the y-ray 
pulse are consistent with both interpretations: they could be the 
result of density bumps swept up by the blastwave (vNiR/ pt < 
v c ), or other shells emitted soon after the first one and catch- 
ing up with the shock front 80 and 140 s after its deceleration 
(Sect. |5J}- 

During the late rebrightening, we know that the spectrum is 
likely described by a broken power law with /3 X = 1.5 + 0.2 
and yS opt = 1 .0 + 0.2. If we interpret it as the cooling break in the 
slow cooling regime, the electron power-law distribution index is 
p = 3.0 + 0.4. The fact that f3 x was around 1 .0 at the beginning of 
the XRT observations can be explained if the cooling frequency 
crossed the X-ray band during the observations. The change in 
the X-ray decay that would be implied amounts to 1/4 and could 
be still compatible with the pre-break slope, ranging from 1.0 to 
1.4, depending on which model one assumes (Sects. [3~3. 1114.1.21 
14.1.31 ). The optical decay gg = 1.06+0.08 is compatible with that 
expected in the case of an ISM environment: a l £ M = 3(p— 1)/4 = 
1 .5+0.3. Differently, the case of a wind environment is ruled out 
(3 cr): a™ d = (3p - l)/4 = 2.0 ± 0.3. 

Alternatively, if during the late rebrightening the bulk 
Lorentz factor has already decreased to T < l/6>j, where 6jj is the 
jet opening angle, the afterglow should already have experienced 
an achromatic jet break: in this case, both optical and X-ray de- 
cay indices are simply equal to p. Notably, this is compatible 
with the measured values during the decay of the late rebright- 
ening and this led Panaitescu] d2007l) to favour the jet interpre- 
tation for this burst. However, we believe that the late rebright- 
ening is more likely to be due to an energisation of the blast- 
wave shock front that strongly affected the measured power-law 
slope, similarly to wha t was inferred in the case of GRB 050724 
dMalesani et al.Ll2007l) . 

The interpretation of the final steep decay following the flare 
as the post jet-break decay seems unlikely: while the late X- 
ray flare might be still interpreted as a shallow-steep transition 
(Sect. [XXXb , in the optical the identification as a flare cannot be 
questioned. Therefore, the steep decay after the peak simply cor- 
responds to the declining part of the flare and not to a jet break. 

5.1. Deceleration radius of the fireball 

The result of Sect. I4.1.4I suggests that the time origin, moved 
to 19.0 ± 2.4 s before the GRB onset, would correspond to the 
explosion time in the detector frame (hereafter DF), i.e. when 
the shell radii are negligible. Let t exp \ = -19.0 s be the explo- 
sion time measured in the DF. In this context we can derive 
some clues on the fireball evolution produced during the first 
shock corresponding to the y-ray pulse. In this picture the shell 
would expand from f e xpi to t y - 0, when it would start emitting 
y rays. The deceleration time fdec would correspond to the peak 
time of the y-ray pulse, i.e. 39.0 + 0.8 s. First the shell accel- 
erates until it reaches the coasting radius (typically ~ 10 13 cm), 
while the bulk Lorentz factor increases linearly with radius un- 
til it reaches the maximum value of To. At this stage, the inter- 
nal energy of the fireball has been converted into bulk kinetic 
energy. After that, the shell expands with constant T = Fo un- 
til it is decelerated by the surrounding medium. If the time it 
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takes the shell to reach the coasting phase is negligible with 
respect to the time it takes to begin to decelerate, it follows 
that during most of the time from f e xpi to t y the shell was mov- 
ing with Fq. In this case the shell would start decelerating at 
the radius 2c(t Y - ? exp i)r 2 /(l + z ) = 1.1 x 10 16 <T /100) 2 (1 + 
z) _1 cm. The deceleration process would culminate at the de- 
celeration radius rd ec = (3 E^IAn nmpC 2 !^) 1 ' - (5.4 X 
10 16 cm) (£ iso 52 /«) 1/3 (r /100)- 2/3 at the deceleration time (DF) 
fdec = 94(£i m .W») 1/3 (rn/100r 8/3 (l + z) s in the th i n shel l 
case ploughing into uniform ISM dRees & Meszarol, Il992[) . 
We note that from the beginning to the peak of the decelera- 
tion, i.e. from t y to fdec, the bulk Lorentz factor decreases ap- 
proximately from To to Fo/2. From this we can infer the dis- 
tance travelled by the fireball in the same time interval: this 
must be between 2c(fd ec - t 7 ) (To/2) 2 /(l + z) and 2c(fd ec - 
t y ) F 2 /(l + Z ), i.e. between 0.6 x 10 16 (r /100) 2 (1 + z)~ 1 cm and 
2.3 x 10 16 (F /100) 2 (1 + z)" 1 cm. This is comparable with the 
distance already travelled at the beginning of the deceleration 
and consistent with the numbers reported above. The unknown 
redshift z unfortunately makes it too tentative to push these esti- 
mates any further. 

6. Conclusions 

GRB 070311 is a FRED-shaped burst followed by early 
NIR/optical pulses detected from ~ 100 to ~ 200 s and possi- 
bly accompanied by a simultaneous faint y-ray tail, with sub- 
sequent variability detected at different NIR/optical filters up to 
~ 10 3 s since the GRB onset. Another remarkable property ex- 
hibited by GRB 0703 1 1 is the late R and X-ray rebrightening ob- 
served around few 10 5 s after the burst, with the X-ray peaking 
earlier than /v-filter photons. When we fit it with the combination 
of an underlying power law plus a pulse, the X-ray (R) fluence of 
the pulse alone is comparable (10 times larger) with that of the 
early pulse, while the peak intensity is about one order of magni- 
tude larger than that of the underlying power law. Interestingly, 
if we refer the times to 19.0 ± 2.4 s prior to the GRB onset, it 
turns out that the peak time as well as the rise and decay times 
of the late pulse are compatible with the corresponding times of 
the y-ray pulse, rescaled by the same factor: / s = 5700 + 700. 
We interpreted this in the external shock scenario, where the y- 
ray prompt emission would correspond to the deceleration of the 
blastwave sweeping up the surrounding medium with uniform 
density, while the late R and X-ray rebrightening would be pro- 
duced by the refreshed shock of another shell emitted after the 
first and impacting the blastwave when this has a Lorentz factor 
F2, so that: T\/75 < T2 < T\. In this context, the time offset 
of 19.0 + 2.4 s before the GRB would correspond to the explo- 
sion time in the detector rest frame, while the GRB onset would 
mark the beginning of the y-ray emission due to the decelera- 
tion of the fireball. From the explosion, onset and peak times 
we infer consistent estimates of the deceleration radius, a few 
xl0 16 (F /100) 2 (l +zy l cm. 

The interpretation of the late rebrightening as the result of 
a density bump in the surrounding medium would explain nat- 
urally both the prompt and the late hump with a single shell. 
However, the possible presence of the cooling break close to the 
X-ray band around the peak time, combined with the remarkable 
flux enhancement observed, AF/F ss 10, makes this scenario 
less appealing. 

The occurrence of the early NIR/optical flares at 80 and 140 s 
after the peak of the y-ray emission is consistent with both sce- 
narios: either due to density enhancements of the matter encoun- 



tered by the blastwave or explained by further shells catching up 
with it. 

Finally, we note that the early (late) NIR/optical pulses do 
not peak contemporaneously with the corresponding y-ray (X- 
ray) pulses, but are delayed by a factor of ~ 3 (~ 2) in time. A 
clear interpretation of this effect is lacking. 
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